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ABSTRACT 

C 17 J = 2 — > 1 observations have been carried out towards the pre-stellar core 
L1689B. By comparing the relative strengths of the hyperfine components of this line, 
the emission is shown to be optically thin. This allows accurate CO column densities 
to be determined and, for reference, this calculation is described in detail. The hydro- 
gen column densities that these measurements imply are substantially smaller than 
those calculated from SCUBA dust emission data. Furthermore, the C 17 J = 2 — > 1 
column densities are approximately constant across L1689B whereas the SCUBA col- 
umn densities are peaked towards the centre. The most likely explanation is that CO 
is depleted from the central regions of L1689B. Simple models of pre-stellar cores with 
an inner depleted region are compared with the results. This enables the magnitude of 
the CO depletion to be quantified and also allows the spatial extent of the freeze-out 
to be firmly established. We estimate that within about 5000 AU of the centre of 
L1689B, over 90% of the CO has frozen onto grains. This level of depletion can only 
be achieved after a duration that is at least comparable to the free-fall timescale. 

Key words: radiative transfer - ISM: globules - ISM: individual: L1689B - stars: 
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1 INTRODUCTION 

Molecular line profiles from pre-stellar and protostellar ob- 
jects potentially offer the best opportunity to extract dy- 
namical information about the collapse process that leads to 
the formation of stars. However, it is becoming clear that the 
interp retation of these li n e pro files can be fraught with diffi- 
culty. Rawlings & Yates (2001) used a self-consistent chem- 
ical and dynamical model of collapsing star-forming cores 
to explore the effects of abundance variations. They showed 
that the line profiles can be very sensitive to the assumed 
values of the free parameters in the chemical models. The 
depletion of molecular species due to freeze-out can have a 
profound influence on the line profiles. Accurate abundances 
are required deep within these cores in order that any freeze- 
out is characterised properly. 

The widely used tracer of molecular hydrogen, CO, is 
so abundant that, in the cold dark clouds in which stars are 
forming, it has a large optical depth and thus cannot trace 
the densest material. Very rare CO isotopomers are there- 
fore used and in order of decreasing abundance these are 
13 CO, C ls O, C 17 0, 13 C ls O and 13 C 17 0, recently discovered 



these conflicting demands but the analysis can be compli- 
cated by the lack of an optical depth estimate. In contrast, 
C 17 is slightly less abundant than C ls O and has a complex 
hyperfine structure revealed in several distinct line compo- 
nents. For a well detected line the hyperfine structure can 
be used to identify whether optical depth effects are present 
and hence abundances can confidently be calculated. 

Of course, in order to calculate hydrogen column densi- 
ties, abun dance ratios b e tween the different isotopomers are 
required. Bensch et al. (2001) demonstrated how by using 



by Bensch et al. (2001) toward the p Ophiuchi molecular has frozen onto the surfaces of grains. Gibb & Little 



several different isotopes of CO, the optical depths could 
be cross-checked and the abundance ratios measured. Ef- 
fects such as isotope-selective photo-dissociation and chemi- 
cal fractionation can be important in translucent clouds but 
are not thought to be s ignificant in the co ld dense environ- 
ments considered here ( Bensch et al. 2001 ). 

Dust emission measurements can be used to calculate 
dust column densities and using a dust to gas ratio, a second 
estimate for the column density of hydrogen can be made. 
If this exceeds the hydrogen column density as traced by 
CO then it is possible that the CO is depleted from the gas- 
phase. The most probable reason for this is that the CO 

~ fl!99^) 



cloud. Measuring abundances is not straightforward however 
because even some of the rare isotopes are mildly optically 
thick and the very rare isotopes can require impractically 
large integration times. C ls O is often selected as a result of 



find CO abundances reduced by a factor of at least 10 to- 
wards the HH24-26 molecular cloud core by this method. 
They also carefully consider alternative explanations that 
do not require abundance differences such as optical depth 
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and b eam filling effects a nd show them to be unlikely. Re- 
cently Bergin et al. (2002) have found that not only CO but 
N2H + is depleted in the Bok globule B68. This raises the 
possibility that in some cores there will be very few avail- 
able infall trac ers. Two very r e cent papers have invesitgated 
CO depletion. Bacmann et al. ( 2002| ) have investigated seven 
pre-stellar cores and find C O to be underab u ndant by factors 
of 4-15 amongst the cores. j0rgensen et al. ( 2002) modelled 
18 pre-stellar cores and found that the class sources in 
their sample are depleted in CO by an order of magnitude. 
L1689B is a pre-st ellar core located in Ophiuchus. 



Gregersen & Evans ( 2000) have observed L1689B in the lines 
of HCO+ J = 3 -> 2 and H 13 CO + J = 3 -> 2. Their mul- 
tipoint data seemed to indicate that the HCO + J = 3 — > 2 
emission was double-peaked and blue-skewed, indicative of 
infall. This is su pported b y the s ingle point CS J = 2 — > 1 
observations by Lee et al] (1999) which sho w a clear blue- 
peak asymmetry. Jessop & Ward-Thompson (2001) observed 
L1689B in C 18 and although they did not have optical 
depth measurements and thus column densities for C ls O 
they explored the dust density and temperature parameter 
space and argued that the CO has to be depleted in order 
to be compatible with the continuum data. 

In this paper, observations of the rare C 17 isotope are 
presented and column densities derived and analysed in an 
attempt to quantify any depletion that is taking place. These 
data are discussed in terms of current models for pre-stellar 
cores in general and L1689B in particular. 



2 OBSERVATIONS 

The observations were carried out at the James Clark 
Maxwell Telescope (JCMT), Mauna Kea, Hawaii on the 
night of 2001 August 19. The C 17 J = 2 — > 1 (224.714368 
GHz) rotational transition was observed using the the het- 
erodyne receiver RxA3. The JCMT half-power beam width 
(HPBW) is 19 arcsec at these frequencies. Typical system 
temperatures were 325 K. Five pointings at 20 arcsecond 



(-20,0), (0,0), (20,0), (40,0) were obtained 



offsets lit (-40,0), 
for C 1( J = : 
tion (16:31:46.90, -24:31:55.98). Jessop & Ward-Thompson 



1, centre d at the epoch 1950 posi- 



Evant (2000) us ed a point i ng po 



24:31:45.00) while Lee et al. (199E 



(2001) and Gregersen 
sition of (16:31:46.98, 
used a pointing position of (16:31:43.6,-24:31:40). The data 
were reduced in the standard manner using SPECX and the 
resulting spectra are displayed in Figure ^ 



3 RESULTS 

Due to the spin of the 17 nucleus, C 17 J = 2 -» 1 is com- 
posed of nine hyperfine components. Figure ^ depicts these 
compon ents and th e ir rel ative strengths (from the tabula- 
1998). The hyperfine components are 
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tion of 

shown with very narrow line widths and also with a broad- 
ening of 0.35 km s _1 to illustrate the typical line shape that 
can be expected to be observed in cold quiescent cores. Since 
the stronger components will be affected relatively more by 
non-negligible optical depths, the shape of the blended line 
can be used to verify that the emission is optically thin in 
this line and transition. Figure H depicts an overlay of the 




Velocity (Km/s) 

Figure 2. Fine structure of the C 17 J = 2 — » 1 transition. The 
narrow spikes illustrate the relative velocities and strengths of the 
hyperfine components. The curve is these components smoothed 
with a gaussian profile with a = 0.15 km s _1 . Both the spikes 
and the smoothed curve are normalised to their peak intensities 



data from the central position of L1689B with the expected 
optically thin line shape. We can conclude from the close 
fit that the line does not display any significant saturation 
effects and is thus of low optical depth. A detailed fit, us- 
ing the HFS routine in the CLASS software package gives 
optical depths in the range 0.5 — 1.5. For such low optical 
depths, the actual value obtained from the software should 
be treated with some caution as the deviations from an op- 
tically thin line shape are very small (for example, the fit in 
Fig. is for very low optical depth) and likely to be strongly 
affected by the noise. In our later analysis, we adopt an op- 
tical depth of ~ 0.7 for all positions because this is consis- 
tent with other recent measurements (N.J. Evans II, 2002 
private communication). The results of the fit are given in 
Table |l] which also lists the line centre, line widths and opti- 
cal depths obtained. The net broadening (due to turbulent 
and instrumental broadening), Vbroad — 0.37 km s _1 . This is 



smaller than the C O li ne width than measured by Jessop & 
Ward- Thompson! (^00l|) who find Vfwhm ^ 0.6-0.8 km s _1 



for their singly peaked and gaussian line profiles (though as 
discussed below, this line is actually optically thick). 

3.1 CO column densities 

If it is assumed that the optical depth is small, local thermo- 
dynamic equilibrium holds, and the excitation temperature 
Tex >> Tomb then the total column density of a species is 
given by 



3k Q(T CX ) / E u 
Ntot = ex P ' 



8n 3 u h 2 S 



(feTex) 



r mb dv: 



(i) 



where J T m \>dV ~ T^AV/r/B, the integrated line intensity. 
If the integrated line strength is measured in km s _1 , fre- 
quency in GHz, fi in debye (1 debye = (1/3) x 10~ 29 C m) 
and if an SI to CGS conversion factor of 4-7reo is used then 
the column density is 

2 Q(Tox) 



iVtot = 1.67 x 10 14 cm 
E 



x 



exp 



T mb dV. 



(2) 
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Figure 1. C 17 J = 2 -> 1 line profiles from offset positions of (from left to right) -40", -20", 0", 20", 40" 



Offset 


Line centre 


Line width 


T~cl7o 


/ T *A dV 


N cl7o 


Mot 




(km s- 1 ) 


(km s- 1 ) 




K km s- 1 


cm 


_o 

cm 


-40" 


3.577 ± 0.004 


0.393 ±0.012 


0.673 ± 0.240 


1.45 


1.0 X 10 15 


1.3 x 10 22 


-20" 


3.592 ± 0.004 


0.378 ±0.013 


1.028 ± 0.242 


1.43 


9.9 x 10 14 


1.3 x 10 22 


0" 


3.526 ±0.004 


0.370 ±0.012 


1.346 ±0.262 


1.35 


9.3 x 10 14 


1.2 x 10 22 


20" 


3.608 ± 0.004 


0.346 ± 0.013 


0.526 ±0.278 


1.34 


9.3 x 10 14 


1.2 x 10 22 


40" 


3.537 ±0.004 


0.356 ±0.011 


1.462 ±0.273 


1.09 


7.5 x 10 14 


9.9 x 10 21 



Table 1. Optical depths, C O column densities and estimated total column densities for each offset position. 




Velocity (Km/s) 

Figure 3. C 17 J = 2 — ► 1 data and expected line shape if 
the emission is optically thin with a turbulent velocity dispersion 
of 0.15 km s" 1 as in Figure U offset to 3.64 km s -1 . Both the 
data and the expected line shape are normalised to their peak 
intensities 



Th e excitation temperature is taken to be 12K, as inferred 
by Jessop & Ward-Thompson (2001). The rotational energy 

and the partition 



levels are assumed to be thermalized at T c , 
function Q(T CX ) for each species is evaluated by interpolat- 
ing betwe en values obtained from the JPL molecular line 

and these are also listed in 



1998) 



database (Pickett et al 
Table 2. Calculating a hydrogen column density from the 
C 17 emission requires knowledge of three conversion fac- 
t ors: the Kb /CO abundance ratio, which we take to be 3700 
( Lacy et al] _ 1994); the oxygen iso tope ratios le O/ 18 0, taken 
to be 560 (jwilson fc Rood||l994|) and 18 0/ 17 0, taken to be 

The final conversion 



3.65 (Ladd et al 



1981) 



|1998| ; [Penziasj 

factor from C 1Y to H 2 is 7.56 x 10° 

Column densities are also calculated using archiv e 
C 18 J = 2 -> 1 data of |jessop fc Ward- Thompson! ( fc00l|) . 



Parameter 


C 17 J = 2- 


+ 1 C 18 OJ = 2^1 


v (GHz) 


224.714 


219.560 


E u (K) 


16.177 


15.8058 


,u (Debye) 


0.11034 


0.11079 


S 


2 


2 


Q(12.0 K) 


4.80 


4.93 



Table 2. Molecular data for C 17 and C ls O. 

This is possible since the C lr O J = 2 — > 1 line is optically 
thin and so the o ptical depth of th e C 18 J = 2 — > 1 line 



can be calculated. Ladd et al 



(1998) show that simply using 



the ratio of the peaks in these lines to calculate the optical 

depth is unsuitable. This is because the closely spaced fine 

structure of the C 17 J = 2 — > 1 line means that the peak 

in this transition is very sensitive to the turbulent velocity. 

The integrated intensities are much more robust and t he ra - 

tio of these are compared with figure 1 of Ladd et al. (1998) 

rm — 



to infer that the optical depth in the C OJ = 2— >1 
line is tciso ^ 2. When the optical depth is moderate as 
in this case, Eqn (|l]) multiplied by a correction factor of 
r/(l — exp — r) can be used to calculate column density. 

The column densities derived from the two lines are 
plotted as a function of offset in Figure ^j. The agreement is 
excellent and both datasets show that the column density is 
remarkably flat across the face of the core. 



3.2 Dust column densities 



Evans et al. ( |200l[ ) obtained SCUBA maps of L1689B and 
carried out a detailed radiative transfer analysis of the dust 
emission which enable them to place constraints on the 
temperature and density distribution in the envelope. The 
gas density distribution was well fitted by a Bonner-Ebert 
sphere with a central density of n c = 1 x 10 6 cm -3 and outer 
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Offset (arcsec) 

Figure 4. Column densities across L1689B, measured using the 
C 17 OJ = 2 — > 1 tran sition (crosses), C ls O data of Jessop & 



Ward-T lomusun _£ 3001 ) (circles) and from SCUBA 850 /im radial 
profile fits of [Evans ct ah ( 2001 ) (pluses). The sizes of the symbols 
may be taken as a rough indication of the observational error. 



radius of 3 x 10 AU . The dust opacities used were from Os- 



senkopl &JHemrmg J 1994 V Using the azimuthally averaged 



results of Evans et al. (2001), and their dust opacity, we cal- 
culate dust column densities, based on an assumed gas to 
dust mass ratio of 100 according to the usual formula 



N = S , 85o/[^«;850A t "lH-B85o(T)] 



(3) 



where Ssso = 0.7 Jy is the flux at 850 microns within a 
33" beam, fl is the solid angle subtended by the beam, 
K850(im = 0.018 cm 2 g _1 is the dust opacity, p = 2.2 is the 
mean molecular weight, ran is the mass of hydrogen and 
S85o M m(T) is the Planck function at T = 12 K. The dust 
column densities are plo tted in Figure | l| and are consistent 



with those obtained by Bacmann et aL[ (2000, 2002) (who 



use a different opacity to c a lculat e their column densities) 
and [Ward-Thompson et al.| (|2002|). 



3.3 Depletion of CO 

It is clear from Figure ^ that the two independent tech- 
niques for measuring the column density produce results 
that differ both in magnitude and distribution. As inferred 
by the above discussions, the exact normalisations of the 
curves in Figure ^ are subject to variation depending on the 
adopted values of dust opacity, isotopic ratios and molec- 
ular abundances. The most important aspect of Figure ^ 
is the flatness of the column density plot from CO com- 
pared with the centrally peaked dust measurements. This 
strongly indicates that depletion is occurring because the 
column density should rise for positions that approach the 
dense centre of the object. The most likely reason for this 
depletion is that the CO has frozen to the surface of the 
du st grains deep with in the core. The best fit dust models 
of Evans et al. (2001) have a dust temperature profile that 



drops from ~ 12 K at the edge of this pre-stellar core down 
to a temperature of ~ 7 K at the centre. The combination of 
the higher gas density and the lower dust temperature is the 
likely reason for a substantial freeze-out of CO in the cen- 
tral regions of L1689B. An alternative explanation for these 
data is that the CO lies along the line of sight and is not 
actually part of the L1689B core. This can be ruled out since 



this would require the core to be completely depleted of CO 
in order that the column density remains roughly constant 
across the face of it; 

compared 



Jessop 
C 18 J = 2 



Andre et al 



Ward- Thompson (2001) compared their 
1 data with earlier mm continuum data of 
Using a model for L1689B with power 



( |1996| ) 

law profiles for the density, temperature and CO abundance, 
they explored the parameter space that would be consistent 
with both datasets. They conclude that the CO could be de- 
pleted by up to 95% but they did not have an optical depth 
measurement for their C ls O data. In contrast, our results 
indicate that the ratio of the column densities towards the 
centre is only around a factor of 3 or so. Of course, the CO 
is unlikely to be uniformly depleted and this factor of 3 rep- 
resents the column depletion. The local depletion is likely 
to be much higher. This can be investigated by plotting the 
column densities that result from models where the radial 
density profile, abundance and temperature are all varied. 
To illustrate this in a very simple case, we use a constant 
temperature model and vary the abundance such that in- 
terior to a certain radius -Rf rcGZC the CO is depleted by a 
factor 0.95 due to freeze-out. Instead of using the Bonner- 
Ebert density distributions which are in convenient t o work 
with, we use the parameters of the fit of Evans et al. (2001) 



profile of the form (see, e 


g. Whitworth & Ward- Thompson 


2001; 


Whitworth & Bate] 


2001) 



p(r) 



PoRq 



(4) 



where po = 1.0 x 10 6 cm -3 is the central density and 
Ro = 750 AU is an inner radius within which the density 
is approximately po. For the purposes here, the difference 
between the Bonner-Ebert and Plummer density distribu- 
tions are of little consequence - both approximate p ~ r~ 2 
in the envelope and are roughly constant close to the cen- 
tre. Figure ^ shows the column density that results from 
this simple Plummer sphere model for a range of values of 
^freeze with the top curve showing the underlying density 
distribution. Comparison with Figure ^ shows clearly that 
values of -Rfroczc of about 40" = 5000 AU would be able to 
reproduce the approximately flat CO abundance with a col- 
umn depletion factor of about the correct degree. One could 
of course investigate the parameter space fully and include 
realistic variations in the freeze-out and temperature but we 
defer that to a later paper. 

The close fit of the C lr O J = 2 — > 1 lines by a model 
with only turbulent velocities (i.e. with no infall or outflow) 
indicates that the CO emitting material is static. In contrast, 
the HCO+ J = 3 -» 2 and H 13 CO+ J = 3 -> 2 data of 



Gregersen & Evans (2000) show double peaked line profiles 
indicating both that the HCO + is self-absorbed and that 
the gas is undergoing bulk gas motions. These emissions 
are likely to originate from within the denser regions of the 
core. The flatness of the abundance points across the core 
measured using C 17 J = 2 -> 1 (Figure §) indicates that 
the CO is depleted from the central regions. Thus one can 
picture that L1689B is composed of a dynamical dense core 
which is depleted in CO but in which HCO + is present. This 
is surrounded by a quiescent surrounding region where both 
CO and HCO + are present. 
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of riH quoted above, but we should be aware that the freeze- 
out timescale is a lower limit and would be larger if the 
sticking coefficient were less than unity, or if (as would seem 
likely) the cores have condensed from a less dense state. We 
may therefore conclude that the level of CO depletion may 
provide a sensitive indicator of the age of cores relative to 
their free-fall times and that a C 17 survey of sources at 
various (early) stages of evolution could provide a powerful 
diagnostic of their dynamical status. 



-50 





r (arcsec) 



50 



Figure 5. Column densities from a Plummer-like sphere density 
distribution with a depletion 'hole' that extends from the centre to 
a radius -Rf rcozc . The local depletion within the hole is such that 
only 5% of the CO remains in the gas phase. The upper most 
curve is with no depletion. The curves below it are for values of 



of between 5 and 45 arcseconds. The curves have not been 



convolved with a telescope beam and hence are more centrally 
peaked than the observational data. 



4 CONCLUSIONS 

Rare isotopes of CO have been observed towards the pre- 
stellar core L1689B. By using the hyperfine structure of 
C 17 J = 2 — > 1 the transition is confirmed to be opti- 
cally thin and emitted by quiescent gas. This allows an esti- 
mate of the total column density as traced by this molecule 
to be made. A comparison of this value with that inferred 
by SCUBA dust emission measurements reveal the CO to be 
depleted in the central regions of this object. The magnitude 
and extent of the depletion is estimated by comparison with 
a simple model of a pre-stellar core with an inner depleted 
region. We estimate that within 5000 AU of the centre of 
L1689B, around 90% of the CO has frozen onto grains. The 
dust temperature at this radius is ~ 10 K (from figure 5 of 



Evans et al. 2001). The sublimation temperature of CO is 
~ 20 K so potentially freeze-out could occur in the outer 
regions of the cloud. In practice of course, the timescale for 
freeze-out in the low density outer regions is long and freeze- 
out will occur only in those portions of the cloud with a high 
local density and low dust temperature. 

For the physical parameters of L1689B, the rate of 



freeze-out of CO is given by (Rawlings et al. 1992) 



nco 



4.57 x W 4 d g a 2 T 1/2 CnuScornco cm" 3 s" 



(5) 



where nn and nco are the hydrogen nucleon and CO den- 
sities respectively, mco = 28 amu is the molecular mass of 
CO, d g is the ratio of the number density of grains to CO 
molecules, a is the grain radius, C is a factor which accounts 
for electrostatic effects, Sco = 1 is the assumed sticking co- 



efficient for CO. Using < d g a 2 >: 
T ~ 10 K yields 

nco = 6.0 x 10 _18 nHnco cm~ 3 s~ 



2.2 x 10" 



and 



(6) 



Using a va lue of riH = 2.8 x 10 cm for L1689B (Bacmann 



et al 



2000) we find that in the absence of other CO forma- 



tion and destruction mechanisms, 90% depletion of CO is 
achieved after ~43,400 years. This is approximately half of 
the nominal free-fall timescale (97,400 years) for the value 
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